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Abstract. We present a first attempt to construct an analytic model for a clumped gas and dust disk and apply 
it to the Galactic Centre. The clumps are described as isothermal spheres partially ionized by the external UV 
radiation field. The disk structure formed by the clouds is described as a quasi standard continuous accretion disk 
using adequately averaged parameters of the discrete cloud model. The viscosity in the Circumnuclear Disk is due 
to partially inelastic cloud-cloud collisions. We find two different solutions for the set of equations corresponding 
to two stable cloud regimes: (i) the observed molecular clouds and (ii) much lighter and smaller clouds which 
correspond to the stripped cores of the observed clouds. It is shown that the resulting physical characteristics of the 
heavy clouds and the disk are in very good agreement with all comparable observations at multiple wavelengths. 
A mass accretion rate of M ~ 10 -4 Mq yr -1 for the isolated Circumnuclear Disk is inferred. We propose that the 
Circumnuclear Disk has a much longer lifetime (~ 10 7 yr) than previously assumed. 
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1. Introduction 

The Galactic Centre is surrounded by a large number of 
gas and dust clouds forming a thick disk {Circumnuclear 
Disk CND) up to a radius of ~7 pcfl. This disk was dis- 
covered by Becklin, Gatley, & Werner (1982) interpreting 
their FIR data as a tilted dust ring. It has its minimum 
emission towards the central compact radio source Sgr A* . 
Subsequently, the CND was investigated by several au- 
thors observing the emission of dust, molecules and atoms. 

- Dust: Mezger et al. (1989), Davidson et al. (1992), 
Dent et al. (1993), Telesco et al. (1996) 

- Molecules: Gatley et al. (1986) (H 2 ), Serabyn et al. 
(1986) (CO, CS), Glisten et al. (1987) (HCN) , DePoy 
et al. (1989) (H 2 ), Sutton et al. (1990) (CO), Jackson 
et al. (1993) (HCN), Marr, Wright, & Backer (1993) 
(HCN) 

- Atoms: Lutgen et al. (1986) (Cn), Jackson et al. (1993) 
(Oi) 

They concluded that the CND has a hydrogen mass of a 
few 10 4 Mq. The disk is very clumpy with an estimated 
area filling factor of <I>a ~ 0.1 and a volume filling factor of 
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1 We assume 8.5 kpc for the distance to the Galactic Centre. 



$v ~ 0.01. The clumps have densities of several 10 5 cm 3 , 
radii of ~0.1 pc and gas temperatures >100 K. A typical 
dust clump has A v > 30 m and M H ~30 M . The physical 
parameters of the clouds in the central 2 pc are listed in 
Jackson et al. (1993). They are partially ionized by the 
radiation of the central Hei star cluster (see e.g. Genzel et 
al. 1996). A low density ionized gas with a density of ^10 3 
cm~ 3 (Erickson et al. 1994) and a temperature of a typical 
Hn region of T-7000-8000 K surrounds the clouds. The 
vertical thickness of the CND increases from ~0.5 pc at a 
radius of 2 pc to ^2 pc at 7 pc. The disk rotates with a 
velocity of ^100 kms -1 which corresponds to a Keplerian 
velocity around a central object of several 10 6 M and 
has a velocity dispersion of ~ 30 kms -1 . It is inclined 
by an angle i ~ 20° — 30° relative to the line of sight 
(LOS). Due to its clumped structure, the CND is usually 
assumed to be a rather short-lived transient feature (~ 10 5 
yr) (Glisten et al. 1987). 

At a distance of ~1.7 pc from Sgr A* it has a sharply 
defined inner edge (Marr et al. 1993). The dust and molec- 
ular line emission drops there by an order of magnitude. 
This defines the outer limit of the Central Cavity (CC) 
which contains the Hn region Sgr A West. It was inves- 
tigated in detail by Lo & Claussen, Lacy et al. (1991), 
Roberts & Goss (1993), Lacy (1994). 
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Three dimensional kinematical models were made by 
Davidson et al. (1992) for the dust emission, Marshall 
& Lasenby (1994) for the molecular line emission, and 
Vollmer & Duschl (2000) for the ionized component. They 
conclude that there is a single plane in which the major 
parts of the CND is located. 

As yet, there are only few attempts to construct a 
physical model for the CND. Wardle & Konigl (1990, 
1993) investigated a continuous smooth disk model for the 
CND including the magnetic field. They succeeded in ex- 
plaining the dust polarization observed by Hildebrand et 
al. (1993). On the other hand, Krolik & Begelman (1988) 
constructed a clumpy disk model for AGNs where cloud- 
cloud collisions are responsible for the energy and mo- 
mentum transport and thus for the viscosity. Cloud-cloud 
collisions dissipate orbital energy resulting a net inward 
drift of the clouds. If the angular momentum of one cloud 
is low enough it can be 'captured by the central object, 
providing its accretion fuel' (Krolik & Begelman 1988). 
They mentioned that the clouds observed by Genzel et 
al. (1985) have a column density comparable to the Jeans 
column density. However, they conclude that sclfgravity 
cannot overcome tidal shear. 

Shlosman & Begelman (1987) pointed out that a nec- 
essary condition for the fragmentation of an externally 
heated disk is that the cooling time must be shorter than 
the orbital time. They argued that this is often the case 
for disks in which the temperature is regulated by dust. 
Furthermore, Shlosman et al. (1990) discussed the possi- 
bility of fueling an AGN by the means of a cloudy disk. 

Given the area and volume filling factor of the CND, it 
is clear that a continuous, smooth molecular disk is ruled 
out. A model of a clumpy disk is needed. In this paper 
we present such a model. First we describe the outlines 
of the model disk physics (section 2). We then give a de- 
scription of the viscosity in the disk and discuss the energy 
dissipation mechanism in section 3. We present the model 
for a partially ionized globule (PIG), give the equations 
for the disk, and treat the physical conditions in section 
4. The results and the verification of the assumptions are 
shown in section 5. Section 6 treats the influence of the 
tidal shear on the clouds. We discuss the mass accretion 
rate in section 7. The conclusions are given in section 8. 

2. The outlines of the model 

We assume that during a short accretion event (At ~ 10 6 
yr) an amount of gas of several 10 4 M is driven into the 
Galactic Centre region at distances less than 10 pc. The 
accreted mass is not uniform but has density fluctuations. 
This clumpy medium is exposed to the ambient UV radi- 
ation field due to the population of young O/B stars in 
the Galactic Centre. Low density regions of the infalling 
gas are evaporated rapidly while regions of higher density 
stay molecular and are heated to an equilibrium temper- 
ature during less than an orbital period. In this way gas 
clouds of different masses and different sizes are formed. 
At a given distance from the Galactic Centre only clouds 



with a central density high enough to resist tidal shear can 
survive. In these clouds the thermal pressure is balanced 
by gravitation. They are assumed to be isothermal. Since 
the UV radiation comes mainly from the Galactic Centre 
their radius in this direction is given by the location of 
the ionization front resulting from the incident radiation. 
The radius of the opposite side is given by the pressure of 
the ionized gas which fills the space between the clouds. 
During the infall the clouds have frequent partially inelas- 
tic collisions. These collisions are highly dissipative be- 
cause the cooling time is much shorter than the duration 
of the collision (Krolik & Begelman 1988). These collisions 
can lead to cloud fragmentation, mass exchange, or coales- 
cence depending on the mass ratio and velocity difference 
of the colliding clouds. If the initial infalling gas has a to- 
tal angular moment which is not zero these collisions will 
lead to the formation of a disk structure, corresponding 
to a clumpy accretion disk, i.e. the CND. Within this disk 
structure clouds lose orbital energy during dissipative col- 
lisions. This results in a net inward drift of the clouds, 
i.e. angular momentum is transported due to cloud-cloud 
collisions (Krolik & Begelman 1988, Ozernoy, Fridman, & 
Biermann 1998). 

We adopt two different views to model the CND as a 
clumpy accretion disk as described above. First, the small- 
scale aspect where the gas clouds are treated as isothermal 
sclfgravitating spheres. For each distance to the Galactic 
Centre the cloud structure is calculated with a given cen- 
tral density, temperature, and external pressure due to 
the ionization front. Second, the large-scale aspect where 
the cloud distribution is smoothed out over the whole disk 
resulting in a continuous accretion disk model. Both mod- 
els share the same temperature distribution and the same 
UV radiation field with respect to their distance to the 
Galactic Centre. The connection between the small-scale 
and the large-scale model lies in the link between the cen- 
tral density of the clouds and that of the accretion disk. 
In addition, the disk viscosity due to partially inelastic 
collisions is taken into account in the large-scale model. 

2.1. The gas clouds 

The gas clouds are exposed to the radiation field of the 
central Hei star cluster, which ionizes the illuminated sur- 
face. The gas in the ionization front flows away from the 
clump and fills up the whole volume of the disk providing 
an outer pressure at the cloud surface which is not illumi- 
nated. The external radiation field is also responsible for 
the gas temperature within the cloud, which is assumed 
to be constant. Consequently, the cloud is modeled as an 
isothermal sphere whose boundary is given at the illumi- 
nated side by the location of the ionization front and at 
the shadowed side by the pressure of the ionized inter- 
cloud gas. For a quantitative estimate we can solve the 
virial theorem for the surface pressure P: 

_ c v M cl T QGMl 
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Fig. 1. The electron density corresponding to the bound- 
ary pressure of the cloud as a function of the cloud radius. 
The dotted lines indicate the electron density at a typical 
cloud radius of 0.1 pc. 

where c v is the specific heat, M c \ the cloud mass, T the 
temperature inside the cloud, r c \ its radius, and is a fac- 
tor of order one. The outer pressure can be translated into 
an electron density, taking T=7000K for the intercloud 
medium (Roberts & Goss 1993). We set M cl = 15 M Q 
and T=100K. The deduced electron density is plotted in 
Fig. [TJ. The electron density at a radius of 0.1 pc is of the 
order of 500 cm~ 3 which compares very well with the ob- 
served value (e.g. Erickson et al. 1994) considering the 
approximate character of the model. 

The CND consists of ~1000 gas clouds. Each of these 
clouds has a mass of ~30 M Q (see e.g. Jackson et al. 1993). 

2.2. The accretion disk 

The clouds can have partially inelastic collisions, which 
are responsible for the energy dissipation and the trans- 
port of angular momentum, i.e. the viscosity (Ozernoy et 
al. 1998). The standard viscosity description does not ap- 
ply in this case, because the turbulence in the CND is 
supersonic. Therefore, we derive a new viscosity descrip- 
tion, which is based on Kolmogorov's theory. This viscos- 
ity depends on the disk height and the energy dissipation 
rate. At the smallest scales the energy is radiated away 
by infrared line emission. Thus, the disk is described fol- 
lowing the standard equations for a smooth continuous 
disk (see e.g. Pringle 1981) with the modified viscosity de- 
scription. In this picture the different components (neutral 
and ionized) of the discrete cloud model are completely 
mixed and the disk central density is directly related to 
the clouds' central density at a given distance from the 
Galactic Centre. 

3. The turbulent viscosity 

In a turbulent medium kinetic energy is transferred from 
large scale structures to small scale structure practically 



without losing energy. So there is a constant energy flux 
from large scales to small scales where the energy is finally 
dissipated. Since the velocity dispersion (Av ~ 20 kms -1 ) 
within the disk is more important than the shear, there is 
no preferred transfer direction. The turbulence can thus 
be assumed as isotropic. In this case the similarity theory 
of Kolmogorov applies (see e.g. Landau & Lifschitz 1959). 
The assumption of a universal Kolmogorov equilibrium 
implies that the kinetic energy spectrum of the turbulence 
depends only on the energy dissipation rate per mass unit 
e and the characteristic size of the turbulent eddy I ~ ^ , 
where k is the wave number. The kinetic energy E(k) is 
related to the mean kinetic energy in the following way: 



E(k)dk 



(2) 



where u is the velocity of the medium. Kolmogorov's the- 
ory yields 



E(k) = Ce3jfc-3 



(3) 



where C is a constant of the order of unity. Considering 
a schematic energy spectrum given by E(k) = 0, for k < 
^turb an d for k > kd and by equation (0) for A: tur b < 
k < kd, we can derive expressions for the dissipative scale 
length Id — k^ 1 , the large scale turbulence scale Zturb — 
^turb' an( i the- large scale velocity v^b- 



(- 3 A)> 



(4) 



where v is the large scale viscosity due to turbulence; 



«turb = <U > - e^tur'b 



(5) 



This leads to a relation between the two length scales and 
the turbulent Reynolds number Re = Vt m b • hutb/v 



(6) 



Consequently the turbulent large scale viscosity v is given 

by 



1 i,l 



(7) 



The macroscopic pressure due to the turbulence is given 

by 



P-P- «t 2 urb ~ P- (e- /turb) 3 



(8) 



where p is the overall density in the turbulent medium. We 
recall that, since there is only negligible energy loss during 
the energy transfer from large scales to small scales, the 
energy dissipation rate per unit mass for large structures 
is the same as that for small structures and is determined 
by the actual physical dissipation process in the smallest 
structures (heating and radiation) . If we can identify this 
dissipative process in a given observed system, it should 
be possible to draw conclusions about the large-scale tur- 
bulence. In the next section we give a possible radiative 
dissipation mechanism in order to get an analytic expres- 
sion for e. 
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3.1. Energy dissipation by radiative cooling 

The clouds in the CND are subject to partially inelastic 
collisions where they lose kinetic energy. During the col- 
lision the clouds' gas is heated, radiating away its excess 
energy in the infrared lines of Ol, Cll, and excited H2 
emerging from the hot outer layer of the cloud. During 
a collision a shock front is formed heating the gas in the 
interacting region to ^1000 K. In a detailed model one 
has to account for this enhanced temperature within a 
small region. In order to treat the energy dissipation in 
our large-scale model we are not interested in the details 
of these non-stationary hot layers but in a smooth sta- 
tionary dissipation rate over the whole disk. We therefore 
make the assumption that the time and space averaged 
dissipation rate in these hot layers can be approximated 
by the continuous infrared line emission of the smoothed 
large scale disk. This point will be further discussed in 
section 5.3. 

At temperatures of ^200 K and densities of several 
10 4 cm -3 , which represent the averaged values of our disk 
model, the Ol and Cll line intensities are comparable (see 
e.g. Tielens & Hollenbach 1985). Here, we will only take 
into account the Cll line at 158 /j,m, assuming that it is 
representative for the radiative dissipation of the turbulent 
energy in the disk. 

The cooling function of the Cll line at 158 (im is 
(Spitzer 1978) 



decided to use only a normalized value of the H/H 2 transi- 
tion region. For the self-shielding its value rj^ 2 follows from 
the balance of formation and destruction of H2 (Burton, 
Hollenbach, & Tielens 1990). 



surface 



(11) 



where Iq — 610" 11 s _1 is the unshielded photodissocia- 
tion rate, (3 = 4 10 5 cm -1 is the selfshiclding parameter, 
(R) = 3 10~ 17 cm 3 s _1 is the H2 formation rate coefficient, 
^surface is the hydrogen density at the cloud surface, and 
Go is the FUV radiation field in units of the equivalent 
Habing (1968) flux of 1.6 10~ 3 ergcm _2 s _1 . This gives 
?7 H2 = 7 (Ar atom i c /N tot ). The factor 7 is determined by 
the condition that the self-shielding becomes the domi- 
nant process for n sur f aco /Go > 100 cm~ 3 (Burton et al. 
1990). For the dust absorption the column density is fixed 
by the UV optical depth t\jv ~ 1- This corresponds to a 



constant column density of N. 



dust 



~ 210 21 cm" 2 , result- 



„dust 



N^JN tot . 



mg m rj 

atomic gas is obtained by 



The total volume fraction of 



of 



Acn = 



AE 



AV PDR At 



= 7.910- 27 d c e~ 92/T nH erg cm" 3 s" 1 ,(9) 



where AVpdr is the volume of the photodissociation re- 
gion (PDR), T is the disk temperature, tih the hydrogen 
density in the disk and dn = "c/"h = 1 j s the frac- 
tion of atomic carbon abundance in the disk with respect 
to the solar one ((nc/nn)© — 310~ 4 ). This is valid be- 
cause at the given density (n > n cr ; t ~ 3 10 3 cm -3 ) and 
at a temperature beyond 200 K the level populations are 
simply determined by their statistical weights (i.e. LTE) 
(Wolfire, Tielens, & Hollenbach 1990). 

The radiative cooling takes place in the outer layers 
of the clouds, i.e. in the PDR, which represent only a 
small fraction of the whole cloud volume. In the large- 
scale model the ionized, atomic, and molecular phases are 
mixed. If the total cloud volume is AV, the modified cool- 
ing function is 



A mod 
CII 



AE 
AVAt 



7.910- 2 V~ 92/T «h ergem^s- 1 



,(10) 



where 77 = AU PDR /AU ~ M at omic/M t ot ^ N &tomic / N tot . 
Anatomic is the mass of atomic gas located in the PDR, M to t 
is its total mass, A^ atom ic is the column density of atomic 
gas, and N to t = vn~ x pd is the total column density of a 
cloud of diameter d. 

The column density of the PDR is determined 
by H 2 /CO self-shielding or dust absorption. As the 
C+/C/CO transition region is difficult to determine, we 



This leads to an energy dissipation rate per mass unit 

: A-f/P - 7.9 10- 27 m p y ot e - 92/T p =:£•/>, (12) 

which can be used in the expression for the viscosity and 
the turbulent pressure derived above. 

4. The detailed model 

4.1. Partially ionized globules (PIG) 

The clouds embedded in the Hn region Sgr A West and 
illuminated by the ambient radiation field are assumed to 
be spherically symmetric, isothermal, and selfgravitating. 
We follow the model of Dyson (1968). The incident UV ra- 
diation ionizes the outer rim and produces an ionization 
front. As only photons above 13.6 eV ionize the hydrogen 
atoms in the ionization front, the FUV field can pene- 
trate further into the cloud. There, these photons cause 
the photodissociation of the H 2 molecules. Both the ion- 
ization front and the dissociation front are governed by 
dust absorption and self-shielding effects. In addition, the 
PDR is preceded by a shock front. Thus, the cloud struc- 
ture has three components, the inner molecular core, the 
predominantly atomic PDR and the outer ionization front. 
Here we are only interested in the ionization front, which 
is assumed to be quasi stationary and D-critical where the 
velocity of the shock front is approximately the neutral gas 
sound velocity squared divided by twice the ionized gas 
sound velocity (see e.g. Spitzer 1978). Due to the much 
higher temperature in the ionization front, the electrons 
are torn away from the cloud, giving rise to an effective 
mass loss. Across the front, the following jump conditions 
apply: 



2c,;' 



Pn 



2c? 



(13) 
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Fig. 2. Illustration of the in the text described PIG model. 
The UV radiation comes from the direction of the Galactic 
Centre. The lower boundary is determined by the gas pres- 
sure of the surrounding ionized low density gas. 



where v is the gas velocity and c is the isothermal sound 
speed. The index n stands for the neutral, the index i 
for the ionized gas component which is located directly 
in front of the ionization front. The outward electron flux 
is assumed to be isothermal and stationary because the 
ionization front moves very slowly into the cloud. The jus- 
tification for all these assumptions is given in section 5.3. 
The solution for a spherically symmetric, stationary flux 
is given by 



r ,Cj,i , < 
— = (— )2 exp(- 

r, v 



Ad 



1 9 

n .cf — v 

— = CX P 2 

n, 2cf 



) , (14) 



where v and n are the velocity and the number density at 
the distance r from the cloud centre. The outer limit of 
the cloud is given by the location of the ionization front 
at r = n. Fig. || illustrates this model. 

In the ionization front the equilibrium between ion- 
ization and recombination is governed by the following 
equation: 



<7j J(l — x) = anx 2 , 



(15) 



where x is the degree of hydrogen ionization, n the hy- 
drogen density, J the number of incident UV photons per 
cm 2 and s, a the recombination coefficient, and Ui the ion- 
ization cross section. The ionization cross section is cal- 
culated using the 'on the spot' approximation. It is only 
valid if the recombination time scale is much less than 
the kinematic timescale. This assumption is justified in 
section 5.3. J(r) is given by 



J(r) = J e- {T * (r)+Td(r) '> 



(16) 



where r g (r) and Td(r) are the optical depths of the gas 
and the dust along a radial line from an outer radius Rb 
to an inner radius r. The optical depth of the gas is: 



,(r) 



n<Ji(\ — x)dr — j — x 2 dr . 



The optical depth of the dust is: 

rflk 

Td(r) = & d I ndr 



(17) 



(18) 



where is the UV cross section of the dust per hydrogen 



nucleus. Differentiating equation 4.1 with respect to r and 
substituting J(r) gives 



J e T9<V W 9 = -an 2 e Td dr . 

Integrating this equation from r, 
J(ri) = riid leads to 



Jo 



Td(ri) 



S Td<T) d 



(19) 

to infinity and setting 
(20) 



This is the final equation governing the radiation transfer 
in the ionization front. 

Since the ionization front velocity is assumed to be 
small compared to the local sound velocity of the neu- 
tral gas, the isothermal neutral condensation is in a quasi- 
equilibrium state. Its structure is described by the Lane- 
Emden equation for isothermal spheres: 



a u 
dx 2 



2 du 

x dx 



(21) 



where x = ^- L (47rGp c ) 2 , p c is the cloud's central density 
and G is the gravitation constant. The radial density is 
given by 



p(x) = p c e 



(22) 



The jump conditions over the D-critical ionization front 
lead to the following expression for the gas pressure: 



p n (x) = p(x)cl 



p c e 



-u 2 



2pi(x) = 2p i (? i 



(23) 



Since the density never drops to zero, an isolated isother- 
mal sphere has no boundary. This changes with an ex- 
ternal radiation field. The outer boundary is given by the 
location of the ionization front at the illuminated side and 
by the pressure of the surrounding ionized low density gas 
at the shadowed side. It should be mentioned here that 
clouds with x > 6.5 are gravitationally unstable (Jeans 
unstable). 

4.2. The disk model 

The disk is assumed to be continuous and smooth. We 
follow Pringle (1981) for the disk equations, replacing the 
a- viscosity prescription by the one described in section 3. 
Furthermore, the pressure due to the turbulent viscosity 
is added to the thermal pressure: 



1 

Re 



£ 3 H 3 p 3 



(24) 
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E = Re £-SH-3p- 



p = Z/H, 



p = pGM(R)(-- 



M 
3^ ' 



R 2 + ^ 



P 



pk B T n 



(25) 
(26) 
(27) 

(28) 



The disk parameters are the height H, the surface density 
E, the mass accretion rate M, and the neutral gas temper- 
ature T n ; fee is the Boltzmann constant. Here £ is defined 
by the equation e = £-p. The second term in equation (28) 
represents the turbulent pressure of the clouds p = /tw 2 urb . 
The radial mass distribution M(R) is given explicitly by 
M (R) = M + M a R§ , where M = 3 10 6 M© is the cen- 
tral mass component and Mo = 1.6 10 6 Mq/pc* describes 
the mass distribution of the stellar content. This is close 
to the findings of Eckart & Genzel (1996). In addition we 
assume that the central density of the disk in the contin- 
uous picture is linearly coupled to the central density of 
the clouds in the discrete picture, i.e. p = (3 ■ pc ■ 

Finally we take into account that the clouds are illu- 
minated from the centre. Therefore the ionization front 
exists only on the illuminated side. As described above, 
there is a flux of electrons away from the clouds. These 
electrons fill up the space between the clouds and thus 
build an intercloud medium which is the Sgr A West Hn 
region. Most of its mass is concentrated within the disk 
plane. Its pressure is responsible for the cloud boundary at 
the shadowed side. The electron density of the intercloud 
medium is calculated using Hi — Ci/fl and pt = flE/Hi, 
where fl is the angular velocity of the clouds at a distance 
R from the Galactic Centre. 
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Fig. 3. Crosses: The Hei stars identified by Genzel et al. 
(1996). Greyscale: Neil line emission (Lacy et al. 1991). 
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Fig. 4. The external radiation field and the inner gas tem- 
perature of the clouds as a function of the distance to the 
Galactic Centre. 



4.3. The physical conditions 

The radiation field in the Galactic Centre seems to be 
made by the central cluster of Hei stars (Genzel et al. 
1996). In Fig. | we show the identified Hei stars together 
with the Neil data from Lacy at al. (1991). The clustering 
seen in Fig. pleads us to the assumption that the produced 
radiation field is approximately spherically symmetric and 
follows a i?~ 2 -law. The Lyman continuum (Lyc-)photon 
production rate within 1 pc is A^yc ~ 10 50 s _1 (Mezger, 
Duschl, & Zylka 1996). Since it is assumed that this radi- 
ation is responsible for the heating of the clouds, the tem- 
perature follows a i?~2-law (for direct heating see Puget 
& Boulanger 1985). The absolute values of the tempera- 
ture are in agreement with molecular excitation models 
(Sutton et al. 1990, Jackson et al. 1993). The FUV radia- 
tion field is G ~ 10 5 (Wolfire et al. 1990) at 1 pc which 
is in excellent agreement with the Lyc-photon production 
rate if one assumes a black body temperature of 35 000 



K for this radiation. The number of UV photons per cm 2 
and s and the gas temperature can be seen in Fig. |4| 

The ambient Hn region has a constant temperature of 
T e =7000 K (Roberts & Goss 1993). The deduced isother- 
mal sound speed in the ionized medium is ci = 7.6 10 5 
cms -1 . The recombination coefficient for hydrogen atoms 
is a = 3.3 10~ 13 cm 3 s _1 , the UV cross section per hydro- 
gen nucleus is tJd = 1-11 10~ 21 cm 2 . The Reynolds number 
is Re — 1000, the critical value for the onset of turbulence 
in laboratory experiments. For the factor j3 between the 
disk's central density and the clouds' central density the 
observed volume filling factor is used: (3 = $v = t^j- 

5. Results and verification of the assumptions 

The problem which must be solved consists of a system of 
6 non- linear equations (equations (20), (24)- (28)) together 
with the Lane-Emden equation (21). It is solved using a 
Newton-Raphson iteration. Given the non-linearity of the 
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Fig. 5. The first set of results for the heavy clouds versus 
the distance to the Galactic Centre, (a) Solid: mean den- 
sity; dotted: density at the cloud surface; dashed: central 
density, (b) Electron density at the outer radius. Solid: il- 
luminated side; dashed: shadowed side, (c) Emission mea- 
sure, (d) Atomic gas fraction, (e) Cloud radius. Solid: illu- 
minated side; dashed: shadowed side, (f) Cloud diameter 
at half maximum density. 



equations multiple solutions are possible. We solved the 
equations with the parameters given above starting the 
Newton-Raphson iteration from different initial values in 
order to recover all solutions. These solutions represent a 
stationary clumpy accretion disk which consists of distinct 
clouds. 

It turns out that we can solve the system of equations 
exactly and that there are two different solutions for the 
radiation transfer with the clouds' central density almost 
unchanged. We refer to the two different solutions as the 
heavy cloud and the light cloud solutions. Since the only 
parameter which is not a priori determined is the mass 
accretion rate, we varied this parameter in the range 10 -6 
Moyr -1 < M < lCT 1 M Q yr _1 . Higher values of M give 
higher central densities and larger disk heights H . The 
adopted value of M = 10~ 4 M Q yr _1 is a compromise 
between a high M leading to high central densities in order 
to resist the tidal shear and a low M giving the observed 
small disk height. Both solutions describe stable gas clouds 
in the sense discussed in section 2. 



5.1. The heavy clouds 

Fig. H and ^ show the results for the heavy clouds. Plot 



||(a) shows the number densities of the neutral gas inside 
the cloud. The averaged density is plotted with a solid 
line, the central density with a dashed line, and the density 
at the cloud surface with a dotted line. It is remarkable 
that the central density has an almost constant value for 
distances greater than 3 pc and that there is more than 
a factor of 6 between the surface density and the central 
density. At a distance to the Galactic Centre of 2 pc the 
central density is about ~ 3.5 10 6 cm~ 3 , corresponding 
very well to the value derived from HCN excitation models 
(Jackson et al. 1993), whereas the surface density has a 
value of ~ 6 10 5 cm~ 3 . Since the gas in the PDR is hotter 
and therefore the gas density is expected to be lower than 
derived in this simple isothermal model, the value stated 
above corresponds well to the value derived from the [Oi] 
63 /im/[Cll] 158 ^m ratio and the H 2 1-0/2-1 S(l) ratio 
(Burton et al. 1990, Wolfire et al. 1990) and from the low- 
J CO line emission (see e.g. Sutton et al. 1990). Thus, 
the HCN line emission originates in the dense cloud core, 
whereas the low-J CO and the H 2 S(l) lines are emitted 
in the atomic/molecular transition region of the PDR. 

Plot |](b) shows the electron density at the cloud sur- 
face as a solid line and the electron density of the in- 
tercloud medium as a dashed line. The latter is in good 
agreement with the value derived from A19 and [Sm] 3 
/im line ratios (Erickson et al. 1994). The electron den- 
sity at the surface lies exactly in the range given by Lo & 
Claussen (1983). We conclude that the radio continuum 
and line emission of high surface brightness at the inner 
edge of the CND originates in the cloud surfaces and the 
low brightness overall emission distribution corresponds 
to the low density intercloud medium. Both features form 
the Hn region Sgr A West. 

Plot ||(c) shows the emission measure EM — 
is in agreement with the value given by 
Beckert et al. (1996) based on the 15 GHz flux of Brown & 
Liszt (1984). The ratio between the column density of the 
atomic gas and the total column density 77 = iV a tomic/-^tat 
(plot ||d) is determined by the dust absorption for dis- 
tances smaller than 2 pc and by H 2 self-shielding for larger 
distances. Estimating the atomic hydrogen column density 
at 2 pc Anatomic ~ 2rjnri gives a value of ~ 410 21 cm~ 2 . 
The observed value (Liszt et al. 1983) is about twice this 
estimate. We observe a significant difference of the cloud 
radius ||(e) between the illuminated and the shadowed 
side. We can interpret them as upper and lower limits 
for the cloud radius. 

The two graphs in plot |^(b) can be interpreted as up- 
per and lower limits for the cloud mass, which lies in a 
range between 10 and 30 M Q . The timescale for the mass 
loss (plot ||c) by the outward electron flux is about 10 8 
yr for a cloud mass of 10 M Q . The clouds are thus in a 
quasi stationary state. An important plot is the one show- 
ing the Lane-Emden parameter (plot ||d) . The dashed line 
corresponds to the shadowed side, the solid line to the illu- 
minated side. Interpreting these two lines again as upper 
and lower limits for the actual x, we find a tendency that 
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Fig. 6. The second set of results for the heavy clouds ver- 
sus the distance to the Galactic Centre, (a) Pressure of 
the ionized gas at the outer rim at the illuminated side, 
(b) Total mass. Solid: calculated with the radius of the il- 
luminated side. Dashed: calculated with the radius of the 
shadowed side, (c) Mass loss rate due to the ionization 
front, (d) Emden Lane parameter. Solid: calculated with 
the radius of the illuminated side. Dashed: calculated with 
the radius of the shadowed side. Horizontal line: limit for 
gravitational collapse x=6.5. (e) Optical depth at 30/im. 
(f) UV optical depth. 



the clouds become gravitationally unstable for distances 
to the Galactic Centre of less than about 2 pc. In order to 
compare our findings with observations, we also plot the 30 
/im optical depth (plot ||e) r 30Mm 2nnK30fj, m where the 
absorption coefficient K3o^m is taken from Draine (1985). 
This can be directly compared with the 30/Ltm measure- 
ments of Telesco et al. (1996), who measured an optical 
depth at 30/im of 10~ 2 at the outer edge of the minispiral. 
This compares very well with the outcome of our model 
(Fig. Us). The UV optical depth across the ionization front 
r uv (plot |6]f) tells us that the ionization front boundary 
is determined by atomic hydrogen self-shielding. 

We now examine the solution for the averaged parame- 
ters of the disk model (Fig. U). The parameters shown are 
the column density E (plot 7a), the disk height H (plot 
0b), the disk's central density (plot 0c), the central pres- 
sure (plot 0d), the velocity dispersion Av of the clouds 
(plot [7)3), the disk flaring angle $ (plot 0f), the radial in- 
ward velocity v r (plot f7g), and the turbulent viscosity v 
(plot fjjh). Since these clouds are only found at distances 
larger than 2 pc, we will only comment on this part of the 



Fig. 7. The disk parameters versus the distance to the 
Galactic Centre, (a) Total gas surface density, (b) Disk 
height, (c) Central gas density, (d) Gas pressure, (e) 
Dispersion velocity of the clouds, (f) Disk flaring angle, 
(g) Accretion velocity, (h) Viscosity. 

graphs. It is remarkable that the disk height fits well with 
the one deduced from observations (Glisten et al. 1987). 
The central density, which is assumed to be proportional 
to the central density of a cloud, stays approximately con- 
stant, which is a reasonable result. The velocity dispersion 
of the clouds lies in a range between 10 and 20 kms -1 , 
compared to 30 kms -1 derived by the same authors. The 
accretion velocity v r is extremely low, which means that 
the clouds take a long time to approach the centre. 

The viscosity v used here is equivalent to an effective 
a as used in the standard prescription v — ac a H. It has 
a nearly constant value of a ~ 0.01. 

The cloud inflow time is given by Krolik & Begelman 
(1988) as 



^inflo 



(Av/v olh ) 2 i coU 



(29) 



where we take the cloud collision rate derived at the 
end of this section, this leads to a mass accretion rate 
of M ~ MtotAinflow Assuming a total disk mass of 
Mtot = 51O 4 M gives M ~ 10~ 4 Mgyr 1 . This value 
is in good agreement with the mass accretion rate used in 
our model. 



5.2. The light clouds 

Figures || and ^ show the same parameters as Figures |^ 
and |6] but for another class of solutions. Here the densities 
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Fig. 8. The first set of results for the light clouds versus 
the distance to the Galactic Centre, (a) Mean density (b) 
Electron density at the outer radius. Solid: illuminated 
side; dashed: shadowed side, (c) Emission measure, (d) 
Atomic gas fraction, (e) Outer radius at the illuminated 
side (f ) Cloud diameter at half maximum density. 



Fig. 9. The second set of results for the light clouds ver- 
sus the distance to the Galactic Centre, (a) Pressure of 
the ionized gas at the outer rim of the illuminated side, 
(b) Total mass calculated with the radius of the illumi- 
nated side, (c) Mass loss rate due to the ionization front, 
(d) Emden Lane parameter calculated with the radius of 
the illuminated side, (e) Optical depth at 30/xm. (f) UV 
optical depth. 



of the cloud centre (plot ||a), at the outer boundary and 
the mean density have the same value than that of the 
heavy clouds. The observed increase of the density with 
increasing distance might not be significant, whereas its 
increase with decreasing distance smaller than 1 pc is im- 
portant because only denser clouds can resist tidal shear at 
smaller distances. The clouds are 100 times smaller (plot 
||c) and 1000 times lighter (plot (9)3) than the heavy ones. 
Thus, the light clouds have the appearance of stripped 
cores of the heavy clouds. 

Their distribution shows a maximum of radius and 
mass at a distance of about 1 pc. The column density 
of atomic hydrogen (plot ^|d) is determined everywhere by 
H 2 /CO self-shielding because of the high density at the 
outer boundary. The electron density at the illuminated 
boundary is 10 times higher than that of the heavy clouds 
(plot |^b), whereas the emission measure is comparable to 
that of the heavy clouds (plot ||c). They evaporate much 
faster (t cva p ~ 10 6 yr at 2 pc) (plot |]c) and are gravita- 
tionally stable (x < 6.5) (plot ||d). 

These clouds might account for the high-J CO lines 
which require a density close to 10 7 cm~ 3 in the 
atomic/molecular transition region (Burton et al. 1990). 



5.3. Verification of the assumptions 

In this section we show that our model assumptions are 
valid for the light and heavy clouds. We only discuss the 
results for the heavy clouds because they contain most of 
the mass of the CND. 

As described above, the outer layer of the cloud has 
three different regions. The ionization front, caused by UV 
photons with energies above 13.6 eV, is located outside, 
followed by the PDR whose origin is the FUV radiation 
field, followed by a shock front due to the higher tempera- 
ture and therefore higher pressure in the PDR. The cloud 
can only be assumed to be stationary if the velocity of the 
shock front u s hock and the velocity of the ionization front 
vif are much smaller than the sound velocity within the 
cloud. In this case they do not affect the clouds' overall 
structure. (Dyson 1968). The velocity of the shock front 
is given by (Bertoldi & Draine 1996) 

« sh ock^ c B (99^*e-^(l + ^)*(l + 0.72g)* , (30) 

o K 

where 8 = Jo/(2an n ri) is the 'Stromgren number', R the 
distance to the Galactic Centre, and tl v is the UV opti- 
cal depth. The velocity of the ionization front can be ob- 
tained by vxe ~ r^M/M where M is the cloud mass and 
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Fig. 10. Heavy clouds properties versus the distance to 
the Galactic Centre: left: the sound velocity in the neutral 
gas (dashed line) together with the velocity of the shock 
front fshock (solid line). Right: the sound velocity in the 
neutral gas (dashed line) together with the velocity of the 
ionization front vif (solid line). 

M ~ AirrjciTii. Fig. |l0| shows these velocities together with 
the sound velocity in the neutral gas. The sound velocity 
is always a factor 10 larger than that of the shock front 
and the ionization front velocity. Thus, the assumption of 
a quasi-stationary equilibrium is justified. In addition, it 
should be noticed that the velocity of the shock front is 
larger than that of the ionization front. 

The assumption of an isothermal ionization front only 
holds if the characteristic cooling length is small against 
the depth of the ionized layer. The characteristic cool- 
ing time in the Hll region outside the cloud is t coo \ = 
6.3 10 11 jm s. With a typical electron density of ~ 10 4 
cm~ 3 we obtain a typical cooling length of L coo i = 
Cit coo \ = 4.8 10 13 cm. This is ~ 10 -4 times the cloud ra- 
dius. If the ionization front has a depth of Ay ~ 0.01 
(Tielens & Hollenbach 1985), its extent is ~ 10 -3 times 
the cloud radius. 

The assumption of an ionization-recombination equi- 
librium is only valid if the timescale of the recombina- 
tion is small against the kinematic timescale. The recom- 
bination timescale at 7000 K is given by Spitzer (1978): 
£rad = 3 10 12 /n o s. With the same electron density used 
above one obtains f ra( j = 3 10 8 s. The kinematic timescale 
is approximately given by t^in ~ r il Q which gives a value 
of t kin ~ 10 11 s. 

Furthermore, we must check if the collisional timescale 
is much smaller than the crossing time for a sound wave 
(thermal timescale). Only in this case does the cloud 
have enough time to reach the quasi-equilibrium state. 
The thermal timescale is given by ithermal = r i/ c n- For 
the collisional timescale we choose the expression v ~ 
^tm-b^coii^ 2 ) (Pringle 1981). In this case it is assumed 
that the collisional timescale is larger than the orbital 
timescale (Shlosman et al. 1990) and the disk can be 
treated in analogy to the Saturn's rings (Goldreich & 
Trcmaine 1978). As the CND has a volume filling fac- 
tor of 0.01 we think that this assumption is justified. 
Nevertheless, it has to be proven by dynamical models, 
which is beyond the scope of this paper. Such a numerical 
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radius (pc) 

Fig. 11. Timescales for the heavy clouds versus the 
distance to the Galactic Centre: Solid line: Collisional 
timescale due to cloud cloud collisions. Dashed line: 
Thermal timescale using the cloud radius of the illumi- 
nated side. Dotted line: Thermal timescale using the cloud 
radius of the shadowed side. 

model, which will be discussed in a following paper, seems 
to agree with this assumption. The different timescales 
are shown in Fig. O. Clearly, the collisional timescale is 
much larger than the thermal one. It is remarkable how 
large this timescale is. In this scenario there is only one 
collision every 10 4 yr in the whole isolated CND. 

We made the assumption that the time and space av- 
eraged dissipation rate in the hot, layers during a cloud- 
cloud collision can be approximated by the Cn line emis- 
sion of the smoothed large scale disk. In the hot partially 
ionized layer which is formed during a cloud-cloud colli- 
sion the gas has a density > 10 5 cm -3 and is heated to a 
temperature of ^1000 K. Under these conditions the Oi, 
Cll, and H2 lines are mainly responsible for the cooling 
of the gas. The total cooling rate in this layer is approx- 
imately A ~ 10~ 17 ergcm _3 s _1 (Tielens & Hollenbach 
1985). The Cn cooling rate of our smooth model disk 
is Aqii = e pc ~ 10~ 19 ergcm _3 s _1 . We now have to 
account for the transient character of the hot layer. Its 
lifetime can be approximated by the thermal timescale 
^thermal of the cloud compared to the timescale between 
two collisions. The transient and continuous cooling rates 
are thus connected in the following way: 



"Oil , ■ V / 

Hhermal 

Fig. [H] shows that icollAthermal ~ 100- This justifys our 
assumption. 

6. Tidal forces 

The clouds in the vicinity of the Galactic Centre are ex- 
posed to strong tidal forces which tend to disrupt the 
clouds, except if self-gravity is effective against tidal shear. 
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Fig. 12. Dimensions of the heavy clouds versus the dis- 
tance to the Galactic Centre. Solid line: the tidal radius 
r c \. Dashed line: the cloud radius at the shadowed side. 
Dotted line: the cloud radius at the illuminated side. 



Fig. 13. Dimensions of the light clouds versus the dis- 
tance to the Galactic Centre. Solid line: the tidal radius 
r c \ of the light clouds. Dotted line: the cloud radius at the 
illuminated side. 



The tidal (Roche) limit for the cloud mass M c \ is given by 
(Mathews & Murray 1987) 

> \M cX rl\f{R)\ , (32) 
5 r c i 5 

where f(R) = GM(R)/R 2 , r c \ the cloud radius, and R 
is the cloud's distance to the Galactic Centre. We have 
calculated the cloud tidal radius r c \ analytically. For the 
heavy clouds it is compared to the cloud radii at the illu- 
minated and the shadowed side (Fig. |l^). The tidal radius 
is smaller than the radius determined by the UV radiation 
field at distances less than 2.5 pc and smaller than the ra- 
dius of the shadowed side at distances less than 3.5 pc. 

We have also calculated the cloud tidal radius r c \ an- 
alytically for the light clouds. It is compared to the cloud 
radii at the illuminated and the shadowed side in Fig. |l3|. 
These very small clouds are stable against tidal shear for 
distances from the Galactic Centre greater than 1.5 pc. 
Comparing the gravitational collapse time tg to the or- 
bital period P rb at a distance of i?c cntrc =lpc from the 
Galactic Centre yields: 

J2_ ^ o^-^-^Y^^J V ° lh A , (33) 
P orb V10 7 cm- 3 / ccntlc Vi50 kms" 1 / 

where n is the gas density of the cloud and u or b is the 
orbital velocity of the cloud. Thus, these self-gravitating 
light clouds will be stretched out along their orbital path 
by an amount comparable to their initial radial extent 
within about one orbital period. 

7. The mass accretion rate 

Glisten et al. (1987) have derived a mass accretion rate 
of the clouds in the CND of M ~ lO^Moyr" 1 . They 
assumed that the total luminosity in the infrared and 
submm-lines comes entirely from the dissipation of the 
turbulent motions. 



We suggest that if the turbulent energy dissipation is 
due to the infrared line emission (for which we have used 
the Cll line as representative), it should not exceed sev- 
eral percent of the whole infrared line emission due to the 
PDRs. The total luminosity in the infrared and submm- 
lines (PDR) is L=2 1O 4 L (Genzel 1989). The emitted lu- 
minosity of the Cn-line L ~ e x Aftot in our model rep- 
resents ~10% of the total luminosity. This is consistent 
with the radiative dissipation rate given by Genzel et al. 
(1985) 

^~"° 3 <£)(^Ki^) 2 <^)W34> 

where rj s is the fraction of kinetic energy converted into 
radiation via shocks, Mcnd is the mass of the CND, and 
At; is its velocity dispersion. We stress here that t co \\ is 
the mean collision time for all clouds in the CND in order 
to recover the total luminosity due to radiative energy 
dissipation. 

8. Conclusions 

We have constructed an analytic model for the CND. It 
mainly consists of ~ 500 heavy clouds moving around 
the Galactic Centre building a disk-like structure. The 
spatial structure of the cloud ensemble is described by 
a smooth continuous disk whose viscosity is due to par- 
tially inelastic cloud-cloud collisions. This work represents 
a first attempt to model a clumped disk analytically in- 
cluding a simple UV radiation transfer. The input param- 
eters (T, Jo, P, Re) are set to observed quantities. At the 
chosen Reynolds number Re turbulence sets in laboratory 
experiments. The real unknown and therefore free param- 
eter is the mass accretion rate M of the disk. The inde- 
pendent resulting parameters (n, n e , EM, ri, H, Av) are in 
excellent agreement with observations at multiple wave- 
lengths. 
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There are two solutions for our set of equations that 
correspond to two stable clump regimes: 

1. the observed heavy molecular clouds with masses of 
~10 M Q and sizes of 0.1 pc, 

2. the stripped cores of the heavy clouds with masses be- 
tween 10~ 5 to 10~ 4 M Q and sizes of ~ 10~ 3 pc. 

It will be of importance to find out whether — in our 
Galactic Center and/or the centers of other galaxies — the 
light cloud population plays a role. 

Within the disk, the number of collisions between 
clouds is very low (n co n ~ 10~ 4 yr _1 for about 500 clouds). 
We infer a mass accretion rate for the isolated CND of 
M ~ lO _4 M yr~ 1 . We conclude that the CND is much 
more stable and has a much longer lifetime (~ 10 7 yr) 
than previously assumed. 

In the present discussion we have focussed on the in- 
vestigation of the two stable regimes of clouds (heavy and 
light). For a more complete model of the CND, however, 
several aspects need to be adressed additionally, including: 

— The influence of tidal forces and internal rotation of 
the clouds on their stability. While tidal forces tend to 
disrupt the clouds, the role of rotation is not so clear, 
since rotation may easily lead to modifications of the 
cloud structure and thus to clouds that are more or 
less stable than the non-rotating ones. The interplay 
between gravity, tidal forces, rotation, and the pres- 
ence of a point-like central mass distribution (i.e., the 
central black hole) may explain why the CND seems 
to be almost empty of material in its inner « 2 pc. 

— The connection between the structure and dynamics of 
the gaseous material in the CND and the stellar com- 
ponent in the immediate vicinity of the central black 
hole. 
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